Some Calculable Contributions to Entanglement Entropy 
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Entanglement entropy appears as a central property of quantum systems in broad areas of physics. 
However, its precise value is often sensitive to unknown microphysics, rendering it incalculable. By 
considering parametric dependence on correlation length, we extract finite, calculable contributions 
to the entanglement entropy for a scalar field between the interior and exterior of a spatial domain 
of arbitrary shape. The leading term is proportional to the area of the dividing boundary; we 
also extract finite subleading contributions for a field defined in the bulk interior of a waveguide in 
3+1 dimensions, including terms proportional to the waveguide's cross-sectional geometry; its area, 
perimeter length, and integrated curvature. We also consider related quantities at criticality and 
suggest a class of systems for which these contributions might be measurable. 



Introduction: The quantum nature of matter is rarely 
evident on macroscopic scales, often due to the decoher- 
ence of excited states toward classical states. However, 
for certain states, such as ground/vacuum states, their 
quantum nature can appear, in principle, on macroscopic 
scales. One of the most dramatic properties of quan- 
tum matter is entanglement and its associated entropy, 
which, if observed on mesoscopic or macroscopic scales, 
would be of broad interest. To define this entropy, con- 
sider a quantum system whose degrees of freedom can be 
divided into two parts in space A, A (see Fig. [ij). The 
geometric or entanglement entropy is defined by the von 
Neumann formula S = —TTj^{pj^ Inp^), where Pa = Tr^p 
is the reduced density matrix of the subsystem A. This 
quantity has appeared in recent investigations in several 
domains including quantum field theory, condensed mat- 
ter physics, quantum computing, and black hole physics. 
It is a measure of one's ignorance of the full system due 
to quantum entanglement between the degrees of freedom 
in the subsystem A and its complement A. 

For d-l-l-dimensional systems with local dynamics the 
entanglement entropy typically obeys an area law S ^ 
A4_i/e'^~^ for d> 2, where Ad_i is the d— 1-dimensional 
area of the boundary dividing the subsystem from its 
complement |T] . Various discussions of this area law have 
been made in the literature, including extensive study 
of bosonic systems in Ref. [5] and fermionic systems in 
Ref . [3] (the latter can involve extra logarithmic factors) . 
Without further refinement the constant of proportional- 
ity is usually ill-defined, as it depends sensitively on an 
ultraviolet cutoff e. By contrast, the entanglement en- 
tropy of 1-|- 1-dimensional systems is well-defined, since 
the e dependence is only logarithmic. For example, the 
entanglement entropy between a pair of half-spaces of 
a 1+ 1-dimensional conformal field theory at correlation 
length ^ was shown to be S* = c/6 In ^/e, where c is 
the central charge of the conformal field theory [4] ; here 
rescaling of e does not alter the coefficient of In^. This 
raises the challenge, to define and calculate unambiguous. 



cutoff-independent contributions to the entanglement en- 
tropy in d > 2 dimensions. 

In this Letter we do just that, albeit in the very spe- 
cial case of free field theory. We place the system in its 
ground state in order to isolate the entanglement entropy. 
Of course real systems are normally at finite temperature, 
which leads to a volume contribution to entropy, but this 
is not our focus. Instead our focus is toward gaining 
insight into novel quantum phenomena at zero tempera- 
ture, such as quantum phase transitions. For a free scalar 
field in d+l-dimensions at finite correlation length ^ (i.e., 
mass fj, = 1/^), we show that in addition to the divergent 
terms, such as 5 ~ Ad-i/e'''~^, there is also a finite area 
law contribution for general smooth geometries 
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for d odd, 
for d even. 
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where 7^ = (-1) V [6 (47r) V((d - l)/2)!]-i for d odd 
and7<i = (-l)'^/2[^2(27r)('^-2)/2(^_j^)|j]-i ford even. For 

a waveguide geometry with specified boundary conditions 
(see Fig. [l] left panel) we define and calculate additional 
power law corrections using heat kernel methods. Those 
methods allow us to express the entropy as an expansion 
in terms of the geometric properties of the waveguide's 
cross-section. We also consider a massless field (^ — >■ 00) 
and define and calculate unambiguous finite terms for the 
interval in a waveguide (see Fig. [11 right panel). 

The area law for general smooth geometries in the mas- 
sive case, as well as the waveguide expansion for both the 
massive and massless cases, extend the results of Ref. [4]. 

Heat Kernel Method: The replica trick is a powerful 
method for computing the entanglement entropy. Since 
the logarithm that appears in the definition of that quan- 
tity is awkward to compute directly, one exploits the iden- 
tity S = -TrA(pAln/9A) = (-|i + l) lnTrp^|„=i and the 
strategy to compute Trp" for integer n and use analytic 



continuation. (A similar analytic continuation is known 
to fail for spin glasses |S], and we suspect that it can fail 
for the entanglement entropy in complicated quantum 
field theories, but it should be safe for the very simple 
field theories considered here.) Consider, for example, 
a field theory in 1 spatial dimension. In that case the 
quantity Trp" is a trace over an n-sheeted Riemann sur- 
face with cut along the subsystem of interest A. If the 
subsystem A is a half-space, then as explained in Ref. pP 
Trp" can be identified with the partition function Z^ on 
a cone of deficit angle 6 — 27r(l — n), and the entropy can 
be recast as 



S= 2Tr 



d6 



l]\nZs 



5=0 



(2) 



which can be calculated analytically. 

Now consider the waveguide geometry shown in Fig. IT] 
left panel. The field lives in the bulk interior of the waveg- 
uide, satisfying boundary conditions on its surface. For 
this geometry we formulate a Euclidean field theory on 
the space CsxA4d~i, where Cs is a 2-dimensional cone of 
radius R (infrared cutoff) and deficit angle S, and Aid-i 
is the {d— l)-dimcnsional cross-section of the waveguide. 
The cone's radius {R -^ oo) corresponds to the physi- 
cal region in space we are tracing over and the angular 
direction is associated with a geometric "temperature" 
(imaginary time) direction in the Euclidean path inte- 
gral. Let Zg be the partition function for a field in its 
ground state defined on this space. For a free scalar field 
of inverse mass ^, the partition function is Gaussian 



InZs = --Indet (-A 



r 



(3) 



where A is the Laplacian satisfying the appropriate 
boundary conditions on Cs x Mj^^i. 

Now let us introduce the heat kernel for the Laplacian 
operator ({t) = tr(e*'^). The trace is defined by im- 
posing Dirichlet or Neumann boundary conditions on the 
waveguide dAAd-i and Dirichlet boundary conditions on 
dCg- This allows us to rewrite the partition function Zg 
and hence the entropy S in terms of the heat kernel: 
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5=0 



(4) 



Since the manifold for the Euclidean field theory for the 
waveguide is a direct product Cs x M.d-l^ the heat kernel 
factorizes as C,{t) — Q{t) C,d-i{t) ■ Thus our problem sim- 
plifies into that of obtaining expansions for two separate 
heat kernels: One for the 2-dimensional cone C,s{t)] and 
the other for the {d — l)-dimensional cross-section of the 
waveguide (^d-i{t)- The heat kernel for the cone has the 
form [6] 



Cs{t) = 
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2tt -5 
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(5) 



FIG. 1: Waveguide geometry in d = 3. Left: Region A is a 
half-space at finite correlation length ^. Right: Region A is 
an interval of length L at criticality. 



where the dots represent terms that either are annihilated 
by the 27r^ -I- 1 1^=0 operator or else vanish in the i? — >■ oo 
limit, and therefore do not contribute to the entropy. We 
thereby obtain 
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1 
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dt 
T 



U-i{t)e-''^^ 



(6) 



for the entanglement entropy for waveguide geometry in d 
spatial dimensions traced over half-space. Evidently the 
entropy is determined by the geometry of the waveguide 
cross-section, through its heat kernel Cd-i(^)- 

Waveguide Cross- Section: The heat kernel for a closed 
domain satisfying either Dirichlet (77 — —1) or Neumann 
(77 = -|-1) boundary conditions in dimensions 0, 1, and 2 
has the following small t expansion ^: 
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Here a is the cross-sectional length of a waveguide in 2- 
dimensions, while A, P, and x are the cross-sectional 
area, perimeter length, and integrated curvature of a 
waveguide in 3-dimensions, respectively. (This expan- 
sion is also of use in computations of the Casimir ef- 
fect between two partitions in a waveguide, see Ref. 1^.) 
The curvature term for an arbitrary piecewise smooth 
2-dimensional cross-section is given by 



X = 



^-^ 24 



IT 



^^ 12tt 



Kilj)d'yj, (8) 



7j 



where ai is the interior angle of any sharp corners and 
k(7j) is the curvature of any smooth pieces. For example, 
X = 1/6 for any smooth shape (such as a circle) and 
X = (»T. — l)/('^ ^ 2)/6 for any n-sided polygonal (so x = 
1/4 for a square). This result differs from Ref. j8] where 
the curvature piece was argued to be proportional to the 
number of corners in an arbitrary shape. 



Regularization - Finite Terms: Direct insertion of the 
heat kernel expansion into eq. ([6]) for the entanglement 
entropy leads to divergences as t — > 0+. These diver- 
gences are associated with the behavior of the theory at 
arbitrarily short distances. As is known, these lead to 
infinities that cannot be renormalized away: logarithmic 
in 1 dimension, linear in 2 dimensions, and quadratic in 
3 dimensions [T]. There are various ways to regulate the 
divergences. For instance, we could impose a hard cutoff 
on the t integral and integrate from t — t^ — (^ to t — go, 
and find terms that only diverge in the e — ^ limit. An- 
other procedure is to use Pauli-Villars regularization by 
subtracting off terms with /i replaced by A and taking 
A large. This is perhaps more appealing as it respects 
the underlying geometry. However, either approach gives 
results containing the cutoff parameters e or A. 

Fortunately, by returning to cq. ([6]) we can identify cut- 
off independent dependence of the entropy on the inverse 
correlation length jjl = 1/^. In general, the leading order 
behavior of the heat kernel as i — )• is 



Q-i{t) 



t{d-l)/2 



(9) 



where a ~ A(i_i/(47r)('^^-'^^/^ is a constant. Inserting this 
into eq. ^ reveals that the integrand has the leading 
order behavior ^ l/i(''+^)/^; giving a divergence of order 
d — 1 as f — > 0+ with respect to a cutoff, say e, defined 
through tc — e^ . (For d — 1 = there is a logarithmic 
divergence.) This singularity can be regulated by taking 
some number of partial derivatives of the entropy with 
respect to the correlation length ^, as that procedure pulls 
down factors of t from the exponential exp(— i/,^^). In 
particular by taking 



Floor 



(10) 



derivatives of S with respect to ^~^ gives a manifestly 
finite integral whose value is independent of any cutoff 
(note k = 1, 1,2 for d = 1,2,3). Hence we define a di- 
mensionless, cutoff-independent quantity through 



-2\fc 



s, ^ i-C') 



d^S 



Using eqs. ^ and ^ and integrating t, we obtain 
1 



(11) 
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for d — 1, 
for d = 2, 



For d = 1 this result is exact [4 . For d > 2 this expansion 
is only valid for a 3> Ci where a is a typical cross-sectional 



length. The third term is topological for d — 3. Note that 
by taking the appropriate number of anti-derivatives it 
is straightforward to isolate cutoff-independent contribu- 
tions to the entropy itself. We have also computed some 
exact results. In d = 2 and for a square cross-section of 
width a in d = 3, Cd-i(0 i^ known exactly. The result 
for S^ is 



acoth(a/^) rj 
24^^ ^ 24' 



for d = 2 



Sf = 
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coth(a/^) + — csch(a/^) 
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- , for d = 3 



where /„^„i = \fr, 



the primed summation means 



{n,™} e Z^/jO, 0}, and Ki is the modified Bessel func- 
tion of the second kind of order 1 . For a 3> ^ we recover 
eq. ( [l2| plus exponentially small corrections. (Note that 
for the square A — a'^,P — ia, x^ 1/4-) 

General Geometries: Although the subleading terms in 
eq. ( 12 1 are specific to a waveguide geometry, the leading 
terms have a meaning for arbitrary geometries. In par- 
ticular, for any boundary in 1 dimension we pick up a 
contribution of 1/12 to S^, as is known [lU]. For closed 
geometries in 2 dimensions, the leading contribution is 
S^ = P/(24^), where P is the perimeter length. For 
closed geometries in 3 dimensions, the leading contribu- 
tion is S^ = A/{A8tt^'^). By integrating up these results, 
we recover the d = 1,2,3 cases in eq. (fTl). Furthermore, 
using the heat kernel in arbitrary dimensions ^ we re- 
cover the general result for arbitrary dimensions. 

This general result differs from estimates made in Sec- 
tion 7 of Ref. [H], where the corresponding term in the 
entropy did not appear. We have checked our result nu- 
merically for the cases of spheres and cylinders, finding 
excellent agreement. In fact our numerics suggests that 
the area term is the only polynomial contribution to S^ 
for large A/^^. We can understand that heuristically, as 
follows: In the regime ^ <C a, where a is a typical length 
scale of curvature of the boundary, the correlations re- 
quired to feel the curvature are exponentially suppressed. 
On the other hand if the boundary contains sharp corners, 
we expect power law corrections to appear. (For related 
discussion, see [H].) We have verified this numerically 
for squares. 

Massless Case - Finite Interval: The previous expan- 
sion requires the field theory to be massive. Let us turn 
now to the critical case (^ — >■ oo). To use our strategy 
to define finite entropy quantities we need a length scale, 
which will now come from considering a finite interval of 
length L, as in Fig. [T] right panel. We can define the 



cutoff-independent quantity: 
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L 



'dL' 



(13) 



The small t heat kernel expansion of eq. ^ is insuf- 
ficient here because we must know the form of Cs (t) not 
only for t < L^, but also for t > L^ where t is large. 
In general the full form of Q{t) is difficult to calculate. 
However, we do not need Csit) for arbitrary S, but only 
the specific limit indicated in eq. El). There are powerful 
tools available for this, as we now explain. The derivative 
of the entanglement entropy can be written in terms of 
an object defined for 2d conformal field theories known 
as the c-function, denoted C. It is related to the inverse 
Laplace transform of ]i^(27r^ -I- 1 ) 0(^)1 5=0 • Convolving 
with the transverse density of states, we have 



Sl^ 



d£C{LV£)pd-i{£). 



(14) 



The c-function C has been studied intensely, see Ref . [13] ■ 
It is known that C(0) = 1/3 and that C is monotonically 
decreasing. Using the heat kernel expansion nh, we can 
inverse Laplace transform to obtain an expansion for the 
density of states pd-i{£). The quantity Sl can then be 
expressed in terms of a few integrals of C, which have 
been computed numerically. We find 



for d = 1, 
for d = 2, 
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for d = 'd. (15) 



Here ki = ^ J^ dx C{x) and k2 = ^ J^ dx x C{x). The 
numerical values are: fci « 0.04 and k2 « 0.01 [11]. For 
d > 2 this expansion is valid for a ^ L, analogous to the 
expansion in (12 1 which was valid for a 3> ^. 

Discussion: (i) We have shown that arbitrary shaped 
domains have an area term, given by eq. (flj), with a cutoff 



independent piece emerging after taking k = Floor[((i -|- 
l)/2] derivatives of S with respect to ^~^. For a waveg- 
uide geometry we used our construction to obtain an 
asymptotic expansion of the entropy for small values of 
the correlation length to cross-section width ratio. For 
arbitrary smooth manifolds the leading order area law 
should be applicable. In contrast to the 1-dimensional 
case, with 5*^ = 1/12, these higher-dimensional entropies 
can be large numerically. It would be interesting to ex- 
tend our results to fields involving alternative dispersion 
relations, fermions, and interacting field theories. 

(ii) Measurement of entanglement entropy of the kind 
discussed above, in the massive, or non-critical, case, re- 
quires changing the correlation length ^ in such a way 
that the microphysics is only weakly affected, and mea- 
suring the corresponding change in entropy AS* (leav- 
ing the cutoff-dependent pieces, such as Ad^i/e'^^^, un- 
affected). Though fluctuations of the vacuum state of 
a relativistic QFT may not be directly measurable J14j . 
we can turn to condensed matter systems. Consider, for 
example, magnetic media. In the absence of an exter- 
nal magnetic field, there is a massless mode; but for an 
externally applied i?-field, (j) acquires an adjustable effec- 
tive mass n — !/£,. In the regime: e^ ^ ^^ ^ A, where e 
is the inter-spin spacing, the area law should be an ade- 
quate description. Another possibility would be to work 
near, but not too near, a quantum phase transition; then 
the correlation length could be varied in a controlled way. 

(iii) A proposal for obtaining the entanglement en- 
tropy experimentally by measuring current fluctuations 
in 1-dimensional electron systems has been presented 
in Ref. [15]. Extending such proposals to 3-dimensions 
would be of great interest. 
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